ESR of coupled spin-1/2 chains: unveiling zig-zag-type interchain correlations 
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By means of electron spin resonance investigations we revealed the crucial role of the interchain 
coupling in the spin dynamics of the spin-1/2 Heisenberg antiferromagnetic (AF) chain material 
copper-pyrazine-dinitrate, Cu(C4H4N2)(N03)2. We found that the dominating interchain interac- 
tion is of a zig-zag type. This interaction gives rise to geometrical frustration effects and strongly 
influences the character of AF ordering. Combining our experimental findings with the results of a 
quasiclassical approach we argue that at low temperatures the system orders in an incommensurate 
spiral state. 



Recently, a large amount of theoretical and experimen- 
tal work has been addressed to low-dimensional (low-D) 
spin systems. Spin-1/2 AF chains can be regarded as one 
of the most important quantum-mechanical model sys- 
tems, exhibiting a rich variety of unusual magnetic prop- 
erties [if. The excitation spectrum of an ideal isotropic 
spin-1 /2 AF chain is gapless and dominated by a so-called 
"spinon continuum" [2J . Since the spin-1/2 AF chain is 
critical, even small perturbations can significantly change 
fundamental properties of the system. One of the most 
prominent examples is a perturbation by an alternating 
g tensor or by the Dzyaloshinskii-Moriya interaction [3j . 
The magnetic excitation spectrum in such a system is 
formed by solitons and their bound states (breathers), 
and the system exhibits a field- induced gap 4]. Inter- 
chain coupling can be regarded as another type of pertur- 
bation, naturally present in all low-D materials. It can be 
considered as a tuning parameter competing with other 
perturbations and determining quantum critical proper- 
ties of spin systems at low temperatures. In addition, 
the interchain coupling can add frustration aspects to the 
problem due to a possible zig-zag-type interaction 043]) 
typical also for spin systems on a triangular lattice. The 
understanding of all these issues is crucial for elucidating 
the low-temperature dynamics of quantum magnets. In 
addition to a comprehensive theoretical analysis, progress 
on the experimental side provides growing opportunities 
to study quantum critical phenomena in low-D spin sys- 
tems. 

The AF chain compound copper pyrazine dinitrate 
(Cu(C 4 H 4 N2)(N03)2, abbreviated as CuPzN) 8] is re- 
garded as the best realization of an isotropic quantum 
spin-1/2 Heisenberg system 0, [l(|. CuPzN has an ex- 
change coupling along the chains of J/fcs = 10.3 K, 



showing no long-range order down to ~ 0.1 K 11]. The 
latter indicates that a perturbation term J', eligible for a 
possible violation of the spin symmetry, is extraordinar- 
ily small. It makes CuPzN an ideal object to investigate 
even tiny perturbation effects in Heisenberg AF spin-1/2 



chains. 

Here, we present results of electron spin resonance 
(ESR) studies of CuPzN, allowing us to probe the quan- 
tum critical behavior of coupled spin chains in the vicin- 
ity of the critical point, separating the Luttinger-liquid 
and the magnetically ordered phases. ESR appears to 
be an effective tool to answer to the question crucial for 
understanding the critical behavior: which type of in- 
terchain interaction is more effective - the transverse or 
zig-zag coupling? The former gives rise to a staggered 
field and a field-induced spin-gap, while for the latter 
the effective field transferred from the adjacent chains is 
cancelled and spin frustration is expected. Analyzing the 
specific ESR properties we studied the influence of possi- 
ble geometrically frustrating interactions in CuPzN and 
proposed a microscopic picture of the magnetic ordering 
in this compound. 

Plate-like single crystals of CuPzN with typical sizes of 
2x1x0.1 mm 3 were used in our experiments. The samples 
were grown by evaporation of a mixture of pyrazine with 
an aqueous solution of copper nitrate [9] . X-band ESR 
experiments were performed at the Zavoiskii Institute for 
Technical Physics of the Russian Academy of Sciences by 
use of the spectrometer Bruker BER 318. High-field ESR 
experiments were performed at the Dresden High Mag- 
netic Field Laboratory (Hochfeld-Magnetlabor Dresden) 
employing a multiple- frequency ESR spectrometer (simi- 
lar to that described in [12j) equipped with a 16 T super- 
conducting magnet and VDI sources of millimeter-wave 
radiation (product of Virginia Diodes Inc.). 

The crystal structure of CuPzN is schematically shown 
in Fig. 1. It belongs to the orthorhombic space group 
Prima with two formula units per unit cell. The cell 
parameters measured at room temperature [8j show very 
slight contraction upon cooling down to 2 K |13| : a - from 
6.712 A to 6.69166(6) A; b - from 5.142 A to 5.10538(4) 
A; c - from 11.732 A to 11.60022(9) A. There are two 
nonequivalent spin chains running along the a axis. The 
CUN2O2 basal planes for two Cu-sites from neighbor- 
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FIG. 1: (Color online) Schematic structure of CuPzN. The 
canting of the basal planes (dashed lines) for nonequivalent 
Cu sites is shown enlarged for the sake of clarity. The intra- 
chain coupling, J, is along the a axis. The possible interchain 
couplings are either between the equivalent Cu sites along 
b (transverse coupling) or connecting nonequivalent Cu in a 
zig-zag manner as shown by the wavy lines. 
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FIG. 2: (Color online) Angular dependence of the (/-factor 
with magnetic field applied in the be plane at 4.2 K (circles). 
The solid black line is the result of the superposition of two 
angular dependences (shown by the dashed and dash-dotted 
lines) originated from two inequivalent Cu sites. 



ing chains are canted by ±5.2° (measured at 4 K) from 
the crystallographic ac plane. Cu-sites from neighbor- 
ing chains are shifted by half a period along a in the 
ac planes, which is not the case for the chains in the 
ab planes. It is reasonable to expect that the quantum 
critical properties will be determined by an interchain 
coupling which is more effective, i.e., either by the trans- 
verse spin-exchange coupling between Cu-sites from the 
unshiftcd chains in the ab planes, or by a zig-zag-type ex- 
change coupling between Cu ions belonging to the chains 
in the ac planes. In the first case, the interchain inter- 
action gives rise to a staggered field and spin-gap, in the 
second case the effective field transferred from adjacent 
chains cancels and spin frustration is expected. 



ESR experiments performed at 9.4 GHz and at temper- 
atures down to 1.4 K revealed a single line with a shape 
very close to Lorentzian; only below 10 K were devia- 
tions from such a shape found. The almost temperature- 
independent (/-factors, g a =2.053, gb = 2.265, and g c 
'2. Ofi.'!. are in g ood agreement with previously published 



ESR data 14J. Figure 2 shows the angular dependence 



of the (/-factor measured at 4.2 K (symbols). Assuming 
tetragonal symmetry, we can describe the experimentally 
obtained angular dependence taking into account contri- 
butions from two inequivalent Cu-sites (shown by the red 
dashed and blue dash-dotted lines). Note that each site 
has the same g\\ and g±, in the corresponding coordi- 
nate system tilted by <j> — ±5.2°) from the basal plane (a 
schematic view along the a axis is shown in the inset of 
Fig. 2). As a result, the (/-tensor for inequivalent sites 
can be described by 



9 = 






2.0543 
±0.0189 





±0.0189 
2.2607 



(1) 



with alternating signs for neighboring chains. 

At temperatures above 10 K, the observed linewidth 
is relatively small (AH ~ 20 — 30 Oe) and shows no 
indication of line broadening due to spin diffusion. On 
the other hand, at lower temperatures, T < J/ks, drastic 
angular-dependent changes in the linewidth, AH , (with 
field applied in the be plane) can be observed (Fig. 3). 
Two maxima located near ±45° appear. (Note, that such 
a behaviour for magnetic fields applied in the ac or ab 
plane does not occur). Interestingly, the ESR linewidth 
exhibits a strong increase at low temperatures when the 
magnetic field is applied in the directions corresponding 
to these maxima (Fig. 4). 

To explain the origin of such dramatic changes let 
us turn to results of high-frequency ESR measurements 
(Fig. 5). The experiments performed at about 55.2, 
110.4, 220.8, and 331.2 GHz allowed us to clearly re- 
solve two ESR lines with g x =2.13(6) and g 2 =2.17(4) 
for magnetic fields applied at 45° with respect to the 
b axis; that corresponds to two inequivalent Cu-sites. 
Note that the g-factor values are consistent to those ex- 
tracted from the X-band ESR data (Fig. 2), for the 
same field orientation where the maximal line broaden- 
ing was observed. A typical angular dependence of the 
high-frequency ESR lineshape is depicted in the inset (a) 
of Fig. 5 which shows ESR spectra obtained at 1.4 K 
and 220.8 GHz with different orientations of the mag- 
netic field in be plane. The ESR absorption is split if 
the magnetic field is applied at 45° with respect to the 
b or c directions (hereafter called I direction), while no 
splitting is observed when the magnetic field is applied 
along the b and c directions. The spectacular tempera- 
ture evolution of the signal is displayed in the inset (b) 
of Fig. 5 where ESR spectra at 1.4 and 42 K for the I 
direction are shown. Two well-resolved ESR lines were 
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FIG. 3: (Color online) (a) Angular dependence of the ESR 
linewidth obtained at 9.4 GHz and at temperatures of 25 K 
(squares), 4.2 K (circles), and 2 K (triangles); the magnetic 
field is in the be plane, b) Angular dependence of the ESR 
linewidth obtained at 9.4 GHz and at 4.2 K for magnetic 
fields in the be plane (circles), ab plane (triangles), and ac 
plane (squares). Lines correspond to fit results (see the text 
for details). 



observed at 1.4 K, which is contrary to the signal at 42 
K. Such a behaviour reflects a crossover between two dif- 
ferent ESR regimes which is driven by the temperature 
dependence of the interchain- hopping rate Uh ■ At higher 
temperatures when the hopping rate exceeds the mis- 
match of the Zeeman frequencies for two nonequivalent 
Cu-sites (|(<7i — g2)\f}H/u)h <C 1; this regime is known as 
fast-fluctuation or line-narrowing regime) only one line is 
observed. Here, /3 = /is/K., where /is is the Bohr mag- 
neton. Upon lowering the temperature, the spin system 
enters the slow-fluctuation regime (| (51— (72)!/? -ff/w/i 3> 1) 
and the two different contributions from two nonequiva- 
lent sites can be resolved. 

Thus, the temperature evolution of the ESR signal, 
shown in the inset (b) of Fig. 5, gives the unambiguous 
evidence of the importance of the zig-zag-type exchange 
interaction (schematically shown by wavy lines in Fig. 1) 
These observations are consistent with results of inelas- 
tic neutron-scattering data @, revealing the absence of 
substantial spin correlations along the b direction. As a 
result, the effective magnetic correlations of CuPzN can 
be described by an anisotropic Hcisenberg model on a tri- 
angular lattice (in the ac plane) with the exchange along 
the a axis (intrachain exchange) equal to J and the weak 
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FIG. 4: (Color online) Temperature dependence of the 
linewidth, AH, measured at 9.4 GHz with magnetic fields ap- 
plied along the b direction (circles) and tilted by 45° towards 
the c axis (squares). Inset: Low-temperature part of the ESR 
linewidth for the 45° orientation; the line corresponds to a fit 
result (see the text for details). 




FIG. 5: (Color online) Upper panel: High-frequency ESR 
spectra with magnetic fields applied along the I direction. 
Lower panel: Frequency-field dependence of the correspond- 
ing ESR absorptions. Inset (a) shows ESR spectra at T = 1.4 
K and 220.8 GHz with magnetic field applied along b, c, and 
Inset (b) shows ESR spectra obtained at 1.4 and 42 K at 
110.4 GHz with magnetic field applied along I. 



interchain zig-zag coupling J± . It was shown recently Q 
that the main features of this model in the limit of small 
J\/J may be captured with its simplest cartoon - the 
zig-zag two-leg spin ladder. 

In the case of a zig-zag spin ladder the spin- 
Hamiltonian can be written as la] 



H = J2[J(SiS i+ i + fif i+1 ) + JxSifr + f i+1 )}, (2) 

i 

where J and J\ are intrachain and interchain exchange 
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interactions, correspondingly, while Si and T, are spin 
operators for two adjacent crystallographic chains. 

The magnitude of the interaction integral Ji may be 
estimated from the analysis of the specific linewidths' 
dependencies at 9.4 GHz (Fig. 3a), which provides one 
more evidence of the triangular interchain exchange. The 
linewidth, AH, may be written as a superposition of two 
contributions 

AH = AH cham + (l/4)( ffl - g 2 ) 2 f3 2 H 2 /io hl (3) 

where the first one is due to the anisotropy of the in- 
trachain exchange, the dipolar coupling, etc. The 
second term describes the inhomogeneous contribution 
which depends on the hopping between the inequivalent 
sites of the adjacent chains. This expression follows from 
the fast-fluctuation limit of the imaginary part of the dy- 
namical susceptibility at the frequency {gx+g*i)l3H/2 [15j 
and has the well-known form of an exchange-narrowed 
linewidth [l6j]. Note also that the first term has a 180° 
symmetry with maxima for magnetic fields aligned along 
6 while the second one, due to the (gi — g-i) 1 factor, has a 
90° symmetry with maxima close to the 45° field orien- 
tation. We assume that the hopping rate (y>h) for T > 10 
K becomes temperature independent and tends to AJi/h 
while for lower T it may be expressed via the dynamic 
scaling relation w t zv [17| . Here, r is the normalized 
temperature, z is the dynamical critical exponent and v 
is the critical exponent of the correlation length. 

To explain the origin of the unusual angular depen- 
dence of the ESR linewidth (Fig. 3) the experimental 
data were described using Eq. (3). For that, gi and 
#2 obtained from the angular dependence at T = 4.2 K 
(Fig. 2) were used. Note, that since at high T the hop- 
ping frequency is temperature independent and equal to 
4Ji/h the fit to the angular dependence of the linewidth 
at T = 25 K gives us J 1 /k B = 0.044 K. Fitting curves 
at this and at lower temperatures, where the exchange- 
narrowing contributions are much more pronounced due 
to the slowing down of L)h with zv ks 1.54, are shown 
in Fig. 3 by lines. Thus, the additional maxima in 
the angular dependence of the ESR linewidth observed 
at low temperatures (T — 4.2 and 2 K, Fig. 3) corre- 
spond to specific directions, where the system approaches 
the boundary separating two types of spin-fluctuation 
regimes {\(g\—g-2)\l3H/uh ~ 1) and the exchange narrow- 
ing process is suppressed. This proves that the critical 
spin dynamic in CuPzN is essentially determined by zig- 
zag-type interchain correlations in the ac plane. 

Let us now turn to the discussion of the possible mag- 
netic ordering of CuPzN, which should strongly be af- 
fected by the zig-zag-type of spin interactions. As men- 
tioned, CuPzN can be described using the effective spin 
Hamiltonian given by Eq. (2). Note, that this model 
is equivalent to a ID spin-1/2 Hcisenberg AF system, 
with nearest-neighbor (NN) and frustrating next-nearest- 
neighbor (NNN) couplings (Ji and J, correspondingly) 



5]. A classical solution for this ID chain system with 
NN and NNN exchange interactions corresponds to an 
incommensurate (IC) spiral structure with pitch angle 8 
Q given by 



cos 9 



-Ji/4J. 



(4) 



For J/k B = 10.3 K and J\fk B = 0.044 K the pitch angle 
is close to 7r/2: 6 = ir/2 + 0.0042. The correlation length 
£, which is given by the relation 



7T/2 ~ 



(5) 



may be estimated as £ ~ 240 in the units of the lattice 
spacing a. Thus, because of the frustrated nature of the 
interchain interactions the ground state in CuPzN cor- 
responds to an IC AF order, while a homogeneous AF 
ordering is impossible at any temperature [3,0]- 

Using the semi-empirical expression for the Neel tem- 
perature [Tij l for a system of coupled spin chains, 



T N = Ach y/]n(\J/T N ) + (l/2)lnln(AJ/rjv), (6) 

with c = 0.233 and A = 2.6 for J x /k B = 0.044 results 
in Tjy = 103 mK. It is worth noting that this value is 
in good agreement with the ordering temperature (107 
mK) found for CuPzN by the use of zero-field muon- 



spin relaxation measurements Uj. This indicates that 
the magnetic ordering observed in llj corresponds to a 
transition to the IC spiral state discussed above. 

In conclusion, quantum critical properties of the quasi- 
1D Hcisenberg AF system CuPzN have been probed by 
means of ESR. The specific angular and temperature de- 
pendencies of the ESR absorption revealed the important 
role of zig-zag-type interchain exchange interactions. Us- 
ing a quasiclassical description and the obtained ESR 
data, CuPzN is predicted to have an IC spiral ordered 
ground state at low temperatures. According to the 
classification discussed in [j| the smallness of the ratio 
J\/J = 0.0043 found in our work makes it possible to 
conclude that CuPzN belongs to an interesting class of 
spin systems which are closely related to the generalized 
Kondo lattice. 
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